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Abstract 

Wc describe the behaviour of the Wilson loops for a novel class of solutions developed 
in [1]. Phase transitions could be found for all source profiles discussed. This study 
has yielded a previously unknown peculiarity not yet seen, where a solution with hi 
hic a profile for which limp_).oo S{p) = 1 can develop a phase transition like behaviour. 
Furthermore, the possibility to create systems with several phase transition has been 
discussed in detail with positive results. 
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1 Introduction 



The AdS/CFT correspondence [2] links a large N d-dimensional conformal field theory to a 
quantum theory in AdSd+i- More specifically, local operators of the conformal field theory are 
related to fields of the dual quantum gravity theory [3] . 

However, AdS/CFT and the more general gauge-gravity theories can be used to study non- 
local operators as well. One of the most interesting of such operators is the Wilson loop [4]. It 
is accessible from the string theory side, which is computationally significant and physically, it 
provides a basis for gauge invariant gluonic operators. 

But even more important here, the potential of a quark-antiquark pair can be related to 
the VEV of a rectangular Wilson loop, whose side are equal to the qq separation and T — )■ oo. 
In the gauge theory the Wilson loop W{C) along a curve C is defined as 

W{C) = ^Tr{P[e'^c^^^dx^]]. (1) 

where is the number of colors, P represents the path ordering of the exponential and 
are the gauge fields. In this paper we will assume the gauge fields to be in the fundamental 
representation, however higher order representations are feasible [6] . On the string side of the 
correspondence we have [5, 7] 

/ W{C)) = [ DFe-^[^l (2) 

^ ' JdF{C) 

where F is used to denote all fields of the string theory and dF their boundary values along C. 

We can approximate (2) using the steepest descent method. Here, the surface spanned by 
the strings ending on C and obeying the Nambu-Goto action Si\ig{F), is minimised. Now, that 
we have approximated the VEV of the Wilson loop and we know that we can relate the qq 
potential to this VEV by ( W{C)) a; e'^^ we see that 

^qq ~ (3) 

Convinced of the importance of Wilson loops, we will compute its properties for backgrounds 
proposed in [1], deepening our understanding of those solutions. These sets of solutions, are very 
interesting as they relate two very fruitful extensions of the original AdS/CFT correspondence. 
One of the extension is referred to as Klebanov-Strassler (KS) Models [16]. Here, a deformation 
is introduced that corresponds to an imbalance in the gauge groups of the field theory side of 
the correspondence. This puts the field theory in the mesonic branch [28, 29]. 

The other extension is usually referred to as wrapped brane models [11, 12]. Here, we obtain 
the field theory through a higher dimensional field theory on which we perform a Kaluza-Klein 
compactification with twisting. 

Our solution are KS like as they share the same IR and UV behaviour, and nontrivially 
relate to wrapped-brane models as they can be constructed from such a model by the addition 
of Nf sources and applying U-Duality, a solution generating technique introduced in [8] . 

2 General Theory 

Here, I will discuss how to compute Wilson loops, as well as the potential energy of the gg-pair 
using the ideas from above for holographic theories. The following is based on [15]. Consider a 
lOd spacetime: 

ds' = -gMid^'f + 9pp{p){dpf + QeMdx'dx^ (4) 
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as well as the Nambu-Goto Action 



Sng — J dadr-^ det[daXi^dj3x''Gi^,^] (5) 

Now using the string gauge r = t and a = x and not allowing the string to explore then internal 
space, we can transform (5) to 

Sng^tJ dxy/P{p{x))+g^{p{x)){d,pY (6) 

Where 

f{p{x)) = guQxx 9\p{x)) = gttQpp (7) 

Now, one can either solve the Eulcr-Lagrange equations of this system or transfer to the Hamil- 
tonian picture and note that the Hamiltonian in the x direction is a constant of motion as the 
Lagrangian is independent of t. Either way, one will derive the following equation of motion 
for the string: 

dp ^.f{p) VP{p)-P{Po) ... 
dx g{p) fipo) ^ ^ 

From here one can calculate that the length L of the static string configuration connecting two 
quarks is 

dx = 2 r^^^Jk^=dp (9) 

Now that we have all the ingredients we need, we can use equation (3) to obtain E{L), the 
potential energy of the quark- antiquark pair as a function of their seperation. This result is 
generally divergent, as we have to assume infinite quark masses in order to make the qq-pair 
static, which in turn allows us to obtain the rectangular Wilson loops needed. The infinite 
potential energy can be renormalized by subtracting the mass of the quarks niq — g{p)dp. 
The renormalized result is 

E = f{po)L + 2 ^^iVm-PiPo) - f{p)) dp - 2 £° g{p) dp (10) 

Now, the backgrounds [1] we will focus on might only be related to wrapped brane models in 
the technical sense discussed in the introduction, but due the fact that wrapped brane models 
are relatively simple to study, we will sketch the construction of our solutions of interest on 
their basis. We start from a generic wrapped brane model, modify it by applying the solution 
generating technique referred to as U-Duahty or rotation and add Nf sources. First, these 
sources are introduced naively, but then will be distributed according to profiles that become 
increasingly complex in order to correct more and more unphysical behaviours that will be 
discussed more detailed later. At every stage we will evaluate equations (9) and (10) to illustrate 
the behaviour of Wilson loops discussed. Afterwards we will construct cases, where the change 
of E(L) to linear behaviour is started off by a phase transition. This is then generalised to 
obtain solutions that produce several phase transitions. 



3 Towards the background 

The wrapped brane models consist of a metric rfs^, a dilaton $ and an RR-three form F3. To 
write down the background, we use SU (2) left-invariant one-forms 

a)i = cos -|- sin sin ^ , 0)2 = — sin -|- cos sin ^ , uj^ — dip + cos 9 d(f , (11) 
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(14) 



and the vielbeins 

E''^ = e^dx\ EP = e^+^dp, E'^ = e'^+^dO , E'^ = e^^'' sinO d^ , (12) 
= ie*+^(a)i + ade), E^ ^ le^+%uj2 - a sin d(p) , E^ ^ ^e^+^uj^ + cos d(fi) . (13) 

Thus, in the Einstein frame, we have 

10 

ds\^Y.{E^)\ 
where we defined 

E'J^-^ = E' NE'^ NE^ N ... NE\ 

h = -27V,e-^-2^ , h = y e-'=-2'^(a2 - 2a6 + 1) , (15) 

3.1 Flavorless Wrapped-D5 

To connect to Section 2, we have to move from the Einstein frame to the string frame. The 
metric can be obtained by multiplying ds'^^ by e 2^ . The metric elements relevant for the Wilson 
loop computation become 

^« = e*('') fc = e*('') Ppp = e*('')+2^(^) (16) 

There are 6 functions in the background (a, b, g, h, k) which depend on the radial coordinate 
p. The background is determined by solving the equations of motion for these functions. A 
system of BPS equations can be derived[9, 17, 18]. These equations show that all background 
functions can be expressed in terms of Q{p) and P{p), where 

Q(p) = 7Ve(2pcoth(2p)-l) (17) 

and P is a solution to the following differential equation 

p»+P'(f;±|:+f;^_4coth(2rt)=o m 

(18) is usually referred to as the master equation. The background functions of interest to us 
are 

^2k_^ 4c&-4^o_ 2sinh(2r)^ 

2 ~ (p2 _ Q2)p/ y^"") 

Interesting solutions can be found semi analytically by numerically connecting the large p 
expansion 



P ^e^^l'' 



c+ + 4p2 - 4p + — + e-^P c_ - — ±p + 



^4g-i6p/3 / 18567 2781 27 



4 " 



/ 18567 2781 27 2 on./sxl 
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with the IR expansion 



hip + 



Ahi 



hi 



P 



IQhi 
^25" 



AN?. 



Shf 3h\ 



(21) 



We solve the master equation numerically, using (21) and its derivative as initial conditions. We 
then obtain the background functions using the relations (19). We then solve the Nambu-Goto 
Action, as outlined in Section 2 to obtain L{pq) and E{pq) given by (9) and (10). These two 
functions can be combined to give us E{L). In this case, this leads to the following plots: 



2 4 6 8 10 12 14 



Po 





Figure 1: We used iV^ = 1, $o = and hi = 2 



As we can see on the left-hand diagram, there is a minimum, non-zero separation length 
beyond which the quarks do not want to be squished together. This gives the field theory a 
characteristic length scale. In other words, the theory has a UV cutoff beyond which we need a 
non-trivial UV completion. This has already been known in [19], which describes the existence 
of an irrelevant dimension 8 operator, which causes the UV of the theory to be not well defined. 
As one can see, for L close to its minimum value, E behaves Coulomb-like, but then grows 
linear for large L, showing that we deal with a confining theory. 



3.2 U-Duality 

So in order to make the connection with Klebanov-Strassler theories, a U-duality as it is de- 
scribed in [8] has to be applied to the system. This procedure is commonly referred to as 
rotation. It is a solution generating technique that schematically unites the backgrounds, by 
taking a solution to our master equation and mapping it to another background where new 
fluxes are turned on. While neither the background functions, nor Q and P change, the metric 
does and new fluxes are turned on. the metric elements relevant to our calculations are 

<s>(p) Hp) , , 

9tt = -^ 9.. = -^ g,, = e''^^^^''^^^Vh (22) 
Vh Vh 

where h = 1 — /t^e^* and n can be chosen to be k = e~'^^°°\ see [20]. The corresponding Wilson 
loop calculations give a very similar picture to the unrotated case. As in the last subsection, 
we still use A^^c = 1 and $o = in the numerical calculations, but /ii = 3 as /ii = 2 gives a 
linear P and unbounded $(p), a solution to which rotation cannot be applied. 



5 



LEE 




L 



Figure 2: The rotated case. The rotation achieved the UV completion needed, which can be 
seen here as limp^,^ooL = 



Evidently different to the last subsection, when we study a solution of P that does not 
remain linear in the IR, but rather grows exponentially L{pq) — > when po ~^ oo. 

The analysis in this subsection has been very similar to [19]. However, in that paper, 
solutions to the master equation (18) have been studied, for which it was assumed that P ^ Q. 
This leads to an IR behaviour of P that is constant to first order and not linear [30] . This has 
led to some interesting effects like phase transitions. For our IR expansion (21) we have not 
been able to find any phase transitions. But in our case, no cusps in the string configuration 
that would break down the steepest descent approximation have been observed either. 



3.3 Flavor branes 

First we add sources in the original, naive way as described in [9, 10], without any kind of 
profile describing the distribution of the Nj branes. This time the form of the metric remains 
unchanged, while the background functions follow different equations. See [20, 21, 22, 9] for 
reference. Q and P are now given by 

P» , H- iV,) [^^^ . - 4coth(2p,] ^ . 

Qip) = \ ^ (2pcoth(2p) - 1) (23) 
We now have for our background functions 

2k ^ P' + Nf 4.i-4#o^ 2sinh(2p)^ 

2 (p2_Q2)(p/ + JV^) ^ ) 

We require new IR asymptotes reflecting the addition of Nf sources. The ones befitting our 
case are [20] 

Pip) = hp + ^ (^-p log p - ^p log(- log p) + (^^Ml^^ ^+Oip' log p) (25) 

Unfortunately, we have adopted an IR singularity caused by the high density stacking of the 
branes near p = [10]. This leads to an unphysical background, which can also be seen through 
the study of the Wilson loops [25] 
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Figure 3: $0 = 0, Nc = l, Nf = 1 and h = 4 

Note the cusp in the function E{L). This is a visual indication of the fact that the way we 
introduced sources causes us problems. 

With sight to the future analysis in this paper, it might be worth mentioning that the addi- 
tion of sources has led our investigation into a regime where we can produce phase transitions 
on the physical (the lower) branch of the E{L) plot: 



LEE 




L 



Figure 4: $0 = 0, iVe = 5, iV/ = 2 and h = 29 



3.4 Sigmoid- like profiles 

As described in [14] , we can avoid this singularity by distributing the sources such that not all 
sources reach p = 0. Of importance is the fact that the density profile for the sources has a 
Maclaurin series with leading term at least of order 0{p^). A very nice and simple profile that 
fulfils the requirements - see [23] - is 

S{p) = tanh^(2p) (26) 

which looks like 

1.0 - 
0.8 - 
0.6 - 
0.4 - 
0.2 - 




Figure 5: S{p) 



7 



Note that its IR expansion is 



S{P) = 16/ - + + 0[p^^] 

as required. Following [14], we now have 

-P' + Q' + 2NfS P' -Q' + 2NfS 



(P" + NfS') + (P' + NfS) 



P-Q 



P + Q 



4 coth(2p) 



(27) 



(28) 



Q = coth(2p) 







coth^(2x) 



-,2k 



P' + NfS 



4$_4(jjQ 



2sinh(2p)^ 



{P^ -Q^){P' + NfS) 



Please note, comparing with the last subsection, it seems that the analysis was simply extended 
by replacing Nf — )■ NfS{p) everywhere. This is not quite true, seen for example at the new 
form the equation governing Q takes. 
We get 




Po 





Figure 6: $o = 0, A^^ = 5, A^/ = 2 and h = 11 



As one can see, E{L) behaves smoothly again. Let us define hic to be the value of hi for 
which P remains linear in the UV. Note that hi ~ hic here, thus P remains linear for quite 
some time, resulting in L{pq) plateauing before approaching 0. 



3.5 Bump-like profiles 



As argued in [1], the UV behaviour is improved by using a profile that decays like e~ . This 
is due to the fact that to preserve the KS like UV asymptotic behaviour, we need a profile that 
decays at least that fast. The presence of an exponentially increasing number of source branes, 
as is the case for a profile where 5 — )■ 1, behaves like the insertion of an irrelevant operator that 
deforms the UV dynamics. We also would not want a profile that decays faster to keep T^"^^^^^^, 
representing the mass density of the sources, positive everywhere. To have — )■ exactly, 
a profile that decays like (sinh(4p) — 4p)~^/^ is needed^. Thus the following two profiles will be 
of importance. See [14, 23] for more details. 

S{p) = tanh^(2p)e^ (29) 

^(P) = ^ uf^l u/3 (30) 
(smh(4p) - 4p)V3 ^ ^ 



thank Jerome Gaillard and Carlos Niiiicz for sharing that information prior to pubhcation. 



8 



which look hke 




0.5 1.0 1.5 2.0 2.5 3.0 



Figure 7: S{p) is blue while S{p) is red 

Note that 5(0) = ^. While this can be quite easily analytically continued to be 5(0) = 0, 
numerical calculations using this profile have to be dealt with carefully. The approach in this 
paper is based on the following idea. As the IR expansion of this profile is 

Sip) = QHAp'-^-fp' + 0[p']) (31) 

we define S{p) = ap^ + hp^ for < p < e. We then solve for a and h by demanding that S{p) 
and S'{p) are continuous at e. For simpler numerics at this stage, let us use S{p) for now, and 
observe that we get a very similar picture to the previous subsections. 



LEE 




L 



Figure 8: $o = 0, A^^ = 1, A^j = 1 and hi = 3.12932 

4 Phase Transition 

An interesting observable phenomenon that can be found is a phase transition. This was first 
observed for SQCD-like field theories in [27]. The interpretation of such a phenomena is similar 
to the phase transition of a boiling Van-der-Waals gas. Several other papers have found such a 
phenomena in similar cases since. See for example [13, 14, 25, 26, 31] 

In order to find them for the theory discussed here, the parameter space of varying hi was 
explored. Please note that not all values of hi are valid. There always exists a critical value 
for hi, For which P grows linearly. For hi > hic, P grows exponentially in the UV and for 
hi < hic, P dies down to a singularity. Note the following figure. 
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Figure 9: The green line represents P for hic, red hic + e and blue h 



Ic 



Our requirement for $ to be bounded only allows exponentially growing solutions like the 
red curve. Thus a shooting procedure was used to find hic, which was used as starting value 
for our exploration into ever increasing values of hi. Unfortunately the search for a phase 
transition was unsuccessful for the three profiles (26), (29) and (30) mentioned above. 

To see how to proceed from here, let us analyse that Van-der-Waals gas analogy from the 
beginning of the section more closely. Such a gas follows an equation of state of the form 



P 



NRT 



V -bN 



(32) 



where R, b, a are constants. Look at the diagrams borrowed from [13] 





Figure 10: The pressure P as a function of the volume V (left panel) and the Gibbs free energy 
G as a function of the pressure P(right panel) for the same isotherm curve. 

Clearly this analogy works by linking L(po) to PiV) as well as E{L) to G{P), such that 
L ■(r^ P, Pq -H- V and E ^ G. On the VdW side, it is clear that a phase transition will be 
observed, if P{V) has a local maximum. This can be achieved by tuning two scales, dictated 
by the constants a and b, representing the interaction between the VdW gas particles and their 
non-zero size respectively. We have so far only explored a 1-dimensional parameter space. To 
increase our chances of success, we have to find another parameter to vary. 



4.1 Introducing 

One thing one can try is to let the source branes not quite reach but only start at some 
point p^. This is easily achieved by multiplying the profile with a step function 0(p — p^) 
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and performing the coordinate transformation p ^ p — p^,. So, for example, our bump profile 
becomes S{p) = Q{p - p*) tanh(2p - 2p,)^e-^(''"^*)/l 

Numerically, one can simply solve the unflavored system between and p* and then connect 
the system by using Punfi(p*) = P{.P*) and P^^fiip*) = P'{p*) as initial conditions for the flavored 
system. As shown in [14] this is justified if S'{p^:) = 5"'(p=k) = which is true in our cases. 

Here, a phase transition can be found regardless of which profile has been used. When we 
studied the parameter space in detail, we found that as we increased hi, L{po) began to flatten 
in the region just below p* until a peak appeared that kept getting more pronounced. See 
Figure 12. This behaviour was so typical that we conjecture that given any acceptable profile 
and p^, > a phase transition can be found for this system. As an example look at the following 
figures. 



L E E 




The position of the phase transition is not random and has a nice physical interpretation. 
Before the phase transition the string explores only the sourcefull region, while afterwards a 
majority of the string is located in the sourceless region p < p*. Thus the phase transition 
appears as more and more of the string enters the sourceless part. 

For a given p^, the value of hi that produces a phase transition is only bound below. This 
also has a nice physical picture associated to it. hi is directly related to c+, the expansion 
parameter of the UV asymptotes of P (20). c+ in turn is related to the scale at which P stops 
being linear and starts growing exponential. This in turn effects the gradient of E{L) and thus 
the "speed" with which the string hits the sourceless region. For a more pictorial description, 
the reader may imagine a light ray getting refracted on a block of glass due to the fact that 
it is slowed down by the medium. The higher hi, the sooner P becomes exponential, and 
the "faster" the string hits the sourceless region. Above a certain threshold we observe phase 
transitions. 



L 




Figure 12: L(po) for the exact same system for hi = hic + 5 at the bottom, till hi = hic + 45 
in increments of 10 and then at the top for hi = hic + 100. 



11 



As the gradient of L{po) grows to oo as po — 0, it is clear that the threshold for producing 
phase transitions, hu also grows to oo as p=K -> 0. This relationship between ^^J^p°^ and hit 
exists as hi represents how much -L(p*) is lifted compared to close by points, and a steeper 
gradient requires more lift. See the figure above. For large po, L is already quite flat, so even 
small values of hi will cause a bump to appear. For small po, L is quite steep, so the initial 
increase in hi only flattens the profile. This is why we have not observed phase transitions just 
varying hi and effectively leaving p^, = 0. 



4.2 Powers of tanh and x 



Of course other variables do also have an effect on the Wilson loop and the associated phase 
transition. 

For example, we can generalise the tanh component in the profiles discussed above from, 
tanh^(2p) to tanh^"(2p) where n > 2 is an integer. This has the conceptual imperfection 
though, that it pushes the flavour correction in the IR expansion of P to higher orders of p and 
thus effectively diminishes the effect of the flavour branes on the system more and more. 



Also X 



being a variable that affects many aspects of this system, also plays a role for 



how easy phase transitions can be found. As an example, here is E{L) for systems equal to 
figure 11. 






4 6 8 10 



7 8 9 10 11 
Figure 13: On the left n = 5, in the middle x = 5 and on the right x = \ 



Please note, that we can also find phase transitions for a; < 1 but higher values of hi have 
to be used. 



4.3 Different profiles 

Again, one could simply take a whole other class of profiles. Very helpful was 



S{p) 



I cosh (4p(jri)+cosh (4pg) 



-De 



-4p 



(cosh (4p)— cosh {4pg)) 2 



^ (cosh (4pqi)— cosh {4pg))^ cosh (4p)- 



1 



2 (cosh (4p)— cosh (4p(j) ) 2 — (cosh (4p)— cosh (4piji))2 
^ (cosh (4p(ji)— cosh (4pi}))-yycosh (4p) — 1 



if p > 4 
if p < pgi 

otherwise 



based on the Flat Measure described in [14]. Here, analytical conditions are found, that the 
profile must fulfil, to produce desirable back reactions. These are then solved in a lengthy 
but logical procedure whose details are mentioned in the paper cited. If one wishes to have a 
decaying profile one could use S'(p) = S{p)e~^''. These profiles have the following form 
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0.5 1.0 1.5 2.0 2.5 

Figure 14: S{p) is blue while S{p) is red 
Both profiles lead to a phase transition. For S{p) we have 




Figure 15: = Nf = 1 and hi = 3.9581 

And for S{p) we get: 




Figure 16: Nc = Nf = 1 and hi = 3.25291 

pqi is related to pq via 

pqi = -arcosh(5 + cosh(4pg)) (33) 

As S decreases, the width of the brane distribution becomes smaller, The mass of the heaviest 
quark decreases and phase transitions are produced. Despite the large amount of control such 
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types of profiles offer, one should be aware, that their realization is not completely clean, in 
the sense that 5" is not continuous everywhere, which can lead to problems as we will see in 
the following sections. 

4.4 A curiosity 

One quite surprising fact that was discovered during the search for phase transitions is that while 
most profiles require hi 3> hic to show transition, profiles of the type S{p) = tanh^(2(p — p*)) 
do not. For hi ^ hic, P{p) will stay linear for quite some time before exhibiting exponential 
behaviour. This can be seen in Figure 9. In such cases, L{pq) will not immediately approach 
0, but first approach a non-zero value a before, and /3 after the transition of P's behaviour. 
Depending on several factors, such as the power of tanh, we can force P > a, leading to a very 
curious phase transition as shown in the following diagram. 



LEE 




Figure 17: Here, Sip) = tanh^(2(p - p*)),iV, = Nf = l,hi = f§f , = 2 

Please note that beside the phase transition the overall shape of E{L) is equivalent to the 
expected, but one would have to plot it approximately between —10'' and 10 to see it. Similarly 
E[p^) peaks at about 10^ before descending rapidly. 

In this section we have described several ways to achieve phase transitions. This can be 
best understood through the variation of two scales, dictated by hi, giving the "speed" of the 
increase of the potential with respect to separation length, and p^,, the scale at which the Nf 
sources kick in. However other techniques as well as influences of the system on the formation 
of phase transitions have been described. In particular, a general characteristic of profiles has 
been identified that is key in producing phase transitions. We need a region of high source 
density transitioning into a region of low density to be able to observe them. This appears 
to be a universal requirement. This obviously leaves the question, if this knowledge can be 
exploited to generalise the profile to describe new and interesting phenomena. 

Also please note that while confinement has not been directly mentioned, it is generally 
understood that the theories discussed are confining. [15] gives several conditions for / and g 
that each guarantee confinement and it is straightforward to check that in all cases at least one 
of them is fulfilled. 

5 Double phase transition 

Something that has been missing so far in the literature, is the possibility to generalise this 
procedure, to find solutions with more than one phase transition. Below we list several types 
of profiles that allow us to generate these. 

In this paper, we will restrict ourselves to profiles producing two phase transitions as proof 
of concept. In general though, these procedures seem to be able to produce an arbitrarily large 
amount of phase transitions if one so wishes. 
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5.1 Double Bump and Tumbling profiles 

The first idea is to add several profiles of the form described in section 3.4 or 3.5. Here, we 
used 

5<(p) _ |tanh'^(2(p - + tanh^(2(p - p^a)) if P > p*2 
ltanh'^(2(p — p*i)) otherwise 

giving us a the following source profile. 



as well as 
S{p) 
looking like 



3 4 5 6 7 



Figure 18: Here, p^i = 1.8 and p^,2 = 5 



tanh^(2(p - p*i))e '^3^*^' 
tanh*^(2(p — p*i))e ^3^^ 



+ tanh^(2(p - p*2))e '"3"*^' if p > p*2 

otherwise 




Figure 19: Here, p=Ki = 1.8 and p=k2 = 5 



This leads to the following plots. For the Tumbling profile we have 
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Figure 20: Here, hi = 47 

And for the Double Bump: 




Figure 21: Here, hi = 47 



Notice, that if the jump is from a region of high source density to a region with lower but 
non-zero source density, the smaller the difference, the harder it is to produce a phase transition, 
independent from the magnitude of the gradient of the jump. This can be seen in figure 20, 
where in the left-hand diagram the second peak is a lot less pronounced than in the right one. 

Also this Tumbling profile in particular can be very interesting for Holographic Technicolor 
theories with Higgsing cascades. This case is discussed in [32]. 

Encouragingly, the central charge of the solutions of all profiles in this section, including 
plateau profiles are monotically increasing and thus show no sign of unphysical behaviour. To 
calculate the central charge, we used the same definition as in [1] 

c = 7^ (34) 

8(apl0g[V/ie2*+2/^+29+fe])3 

This leads to graphs of the following shape: 




12 3 4 

Figure 22: This is the central charge for the Tumbling case. For the other cases, the shape is 
similar, but c increases much faster 



16 



5.2 Plateau profiles 



So far all the phase transitions observed come from the string moving from regions of high 
source density to regions of lower density. It would be interesting to see if this also works the 
other way around. Thus we would like to engineer profiles that quickly rise from to 1, stay 
at 1 for some time before quickly decreasing back to 1. It is a major challenge to find such 
functions and so far none have been found that have a continuous first derivative. These can 
still be used, given that the discontinuity is well away from the regions of interest, such as the 
middle of the plateau, but will lead to an imperfection discussed below. 

The easiest way to obtain profiles as described, is to use a profile that behaves like S{p) 
from section 3.4 and 4.1 and then mirror it around the axis p = c, where S{c) ~ 1 and then 
forcing to be after it reached that value again. As an example we used 

{0 if p > p* + 4 

tanh^(2(p - p, - 4)) if p* + 4 > p > p, + 2 

tanh'^(2(p — p*)) otherwise 

resulting in 




3 4 5 



Figure 23: Here, Nj = = I, = 1.8 and h = 47 
This leads to the following situation. 




0.015 0.03 0.045 7.5 8.0 8.5 9.0 9.5 10.0 



Figure 24: -L(po) and both physical phase transitions 

Please note that E{L) has 3 phase transitions, the two physical ones shown above, and one 
unphysical that is due to the discontinuity of S" at p = 3.8 at which the numerics produce the 
kink seen for L{pq). Note that the background does not contain any singularities per se. It is just 
at the exact point of the discontinuity that the numerics break down. Thus this discontinuity 
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is usually not noticeable and does not produce a third unphysical phase transition. In order to 
highlight this breakdown, L{po) has been evaluated only at discrete points and then interpolated 
inbetween. Furthermore L{po) was evaluated at exactly 3.8 to produce an erronous data point. 
The purpose of this was to highlight the limitations of the numerical calculation, and illustrate 
the dangers of using profiles without a continuous derivative. Also due to the different scales 
and sizes at which the two phase transitions occur, we had to graph them separately to seem 
them clearly. 

Thus we have strong evidence that phase transitions also occur when entering region of 
higher source density. This leads us to believe, that every bump profile actually produces two 
phase transitions per bump, that merge into another through an effect that is analogous to 
the effect that limits the angular resolution of lenses. Both bumps in L{po) or -E'(po) of, for 
example, figure 11 are so close together that they merge into one. 

6 Conclusions 

The basic idea of this paper was to study the details of Wilson loops of a novel class of solutions 
[1], that generalise the KS [16] and baryonic branch [24] solutions, that, by the addition of 
sources, move the QFT to the mesonic branch.. 

This was done constructively, by starting from a relatively plain wrapped branes model and 
slowly adding all the necessary features, keeping track of the properties of the Wilson loops of 
the theory at each stage. 

Then, since a specific profile has not yet been obtained from a microscopic kappa-symmetric 
solution to the brane equations of motions, several profiles that feature desirable qualities have 
been examined. 

Afterwards wc discussed the possibility of creating a phase transition by tuning two scales 
dictated by hi and p^,. Also the effect of other variables of the theory on this phenomenon has 
been described. This construction of phase transitions has been very successful and yielded a 
positive outcome in every case. Furthermore, an interesting and so far unobserved phenomenon 
has been discovered and described. 

Lastly, the possibility to create several phase transitions has been discussed, and several 
features a profile could have to produce these have been examined in detail, yielding a double 
phase transition in each case. In particular, a profile that might become useful for some extended 
holographic technicolor theories has been described. 

My procedures should also be apphcable to the 2+1 d equivalent that has been discussed 
in [26] 

An interesting open question would be a detailed mathematical derivation of the exact 
requirements needed to obtain a phase transition. We may have some results on this in the 
near future. 
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